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We have studied the influence of nearest- neighbor (NN) repulsions on the low frequency phases 
diagram of a quarter-filled Hubbard-Holstein chain. The NN repulsion term induces the apparition 
of two new long range ordered phases (one Akp CDW for positive (7e// = U — ig^ /lo and one 
2kF CDW for negative U^fj) that did not exist in the V^ = phases diagram. These results are 
put into perspective with the newly observed charge ordered phases in organic conductors and an 
interpretation of their origin in terms of electron molecular- vibration coupling is suggested. 



I. INTRODUCTION 

It is well know that low dimensional systems are 
susceptible to structural distortions driven by electron- 
phonons interactions. The most commonly studied 
phonons-driven instability is the metal to insulator 
Peierls transition in one-dimensional (ID) conductors. 
The insulating state is a periodic modulation of bonds 
charge density (BDW) associated to a lattice distortion. 
Such instabilities are driven by the coupling between the 
electrons and the inter-site phonons modes, the interac- 
tion been essentially supported by a modulation of the 
hopping integrals between two nearest neighbor (NN) 
sites. The consequence is an alternation of the bond or- 
ders while the on-site charge remains homogeneous on 
the whole system. 

When a "site" stands for a complex system with inter- 
nal degrees of freedom, there is another important type of 
electron-phonons (e-ph) interaction, namely the one that 
couples the electrons with the internal phonons modes 
of each "site" . This is, for instance, the case in molecu- 
lar crystals where a site represents a whole molecule. In 
such cases the totally-symmetric (Raman active) molec- 
ular vibrational modes caii couple to the system elec- 
tronic structure. HolsteinEI was one of the first to un- 
derstand the importance of such e-ph coupling, showing 
that it may lead to a self trapping of the conducting elec- 
trcms by local molecular deformations. In the 60's Lit- 
tlea even suggested that intra-molecular vibrations could 
be responsible for super-conductivity in the organic con- 
ductors. More recently it was again proposed to be t 
super-conductivity mediator in fuUerides-based system, 

Even-though the electron molecular-vibrations (e-mv) 
has now been excluded as super-conductivity mediator 
in organic conductors, a simple analysis shows that they 
should in anyway be relevant to these systems. In- 
deed, the conductivity-supporting molecules, such as the 
TMTTF or TMTSF -based molecules, the TCNQ- 
based molecules, the M((imji)2-based molecules, etc, 
have a certain number of characteristics in common. 
They are large, planar, conjugated and based on organic 
cycles, all characteristics favorable to a strong coupling 
of the conduction electrons (which, belonging to the tt 



conjugated system of the molecules, are strongly delocal- 
ized on the molecular skeleton) with the molecular to- 
tally symmetric (Ag) vibrational modes. This analysisjs 
fully supported by Raman spectroscopy measurementsQ'El 
which assert both the existence of low frequency vibra- 
tional modes (associated with ring angular deformations) 
and e-mv coupling constants belonging to the intermedi- 
ate regime. 

It is widely admitted that the simplest pertinent model 
for describing the ID organic conductors electronic struc- 
ture is the extended Hubbard (eH) model with NN bi- 
electronic repulsions. The present study seeks there- 
fore at studying the combined effects of the electron 
correlation within the eH model and the e-mv interac- 
tions, in a quarter- filled ID chain, relevant for organic 
ID conductors. The e-mv problem have been largely 
addressed in the case of the one-dimensional half-filled 
chainB. Quarter-filled systems have been treated in sev- 
eral regimes. In the weak coupling regime renormaliza- 
tion group (RG) approachesEl sh«w,that the transition 
line between theiuttinger LiquidEHl3 (LL) phase and the 
Luther-EmeryoE£l (LE) phase (gaped spin channel and 
dominating 2kp charge fluctuations) is displaced toward 
positive gi parameters when the e-mv coupling increases. 
In addition, the Luttinger Liquid parameters are renor- 
malized by the e-mv interactioris,, In the adiabatic and 



small inter-site repulsion regimeE2ltll, small systems diag- 
onalizations exhibit three different phases, one uniform 
phase at small e-mv coupling, associated with LL, one 
2kp charge density wave (CDW) phase for large enough 
e-mv coupling and small values of the on-site electron 
correlation, and one Akp CDW phase for large enough 
e-mv coupling and on-site repulsion p4Vhen inter-site re- 
pulsion is omitted, we have studiedll3 the whole phases 
diagram as a function of both the phonons frequency, the 
e-mv coupling and the on-site correlation strength. We 
have shown that the dependence of the phases diagram to 
the phonons frequency is crucial. Indeed, while for high 
frequencies (corresponding to the highest Ag molecular 
vibrations of the Bechgaard salts) the phases diagram 
is very poor and well reproduced by the weak coupling 
approximation, for low phonons frequencies (correspond- 
ing to the lowest Ag molecular vibrations of the Bech- 



gaard salts), the phases diagram is on the contrary very 
rich. At small e-mv coupling, and in agreement with the 
RG results, we found LL and LE phases with renormal- 
ized parameters. In the intermediate coupling regime, we 
found at surprisingly small values of the on-site repulsion 
(from U/t '^ 2), a metallic phase with dominating Akp 
CDW fluctuations. For large e-mv coupling we found 
polaronic phases where the electrons are self-trapped by 
the molecular deformations, either by pairs (low electron- 
correlation regime) or alone (large electron-correlation 
regime) . 

The NN bi-electronic repulsions are crucial for a re- 
liable description of the ID organic conductors. The 
present paper will therefore study the interplay between 
the e-mv and the NN bi-electronic repulsion within an 
extended Hubbard model. In regards to the previous 
results we will limit us to the low phonons frequencies 
where we can expect significant effects to occur. 

The next section will be devoted to the model descrip- 
tion and the computational details. Section 3 will present 
the results and section 4 will discuss their relevance to the 
organic conductor physics. The last section will conclude. 



II. COMPUTATIONAL DETAILS AND MODEL 
A. Model 

The simplest way to couple dynamically dispersion- 
less molecular vibrations to the electronic structure is 
through local harmonic oscillators and linear e-mv cou- 
pling. We will therefore use an extended-Hubbard- 
Holstein model (eHH) . If U stands for the on-site repul- 
sion, V for the nearest-neighbor coulomb repulsion and 
g for the e-mv coupling constant the eHH model can be 
written as He + Hph + He-mv with 

He = t y^ (4^i.g.c,,o. 4- h.c.) + £/ y^ '^i,T"i,X + ^ y^ "■»"t+i 

Hph = LuJ2iblb, + l/2) 
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Ci^a and rii.o- are the usual creation, annihilation and 
number operators for an electron of spin a located on 
site i {rii = Ui^^ -t-n^^j). bj and bi are the intra-molecular 
phonons creation and annihilation operators and lu the 
phonons frequency. Jhe energy scale is fixed by t = 1. 

As noticed in ref.td, the on-site part of the Hamiltonian 
can be rewritten (apart from constant terms) as 



• The on-site bi-electronic repulsion term is renor- 
malized by the e-mv coupling and the effective in- 
teraction Ueff = U — '2g'^/u) becomes attractive in 
the strong coupling regime. 

• One sees from eq. |l| that the phonons and e-ph 
parts of the Hamiltonian can be rewritten as a dis- 
placed harmonic oscillator. The noticeable point 
is that the displacements are proportional to the 
sites charge, simulating this way the relaxation of 
the molecular geometry as a function of the ionicity 
of the site. 

• The natural basis for the phonons states is therefore 
the eigenstates of the displaced oscillators. Such a 
vibronic representation is not only very physical, 
but also particularly suited for the representation 
of the low energy physics. Despite the fact that 
one would have a complete basis set for each value 
of the sites occupation, the necessity to work in a 
truncated representation lift the problem of over- 
completeness. 

• One should also notice that there is a strong renor- 
malization of the hopping integrals between the ini- 
tial and final vibronic states of two neighboring 
sites, by the Franck-Condon factors. In this rep- 
resentation the Franck-Condon factors correspond 
to the overlaps between the vibronic states asso- 
ciated with ±1 occupation numbers. As physi- 
cally expected, when an electron hopes between two 
sites vibronic ground states, the hopping integral is 
exponentially renormalized by the displacement : 
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On may highlight the following points. 



The direct consequence is 
an increased tendency to electron localization. In- 
deed, vibronic high energy states need to be sum- 
mon up for delocalization processes when the c-mv 
coupling is large. 



• The pertinent e-mv coupling parameter is not g but 
rather g/u;, in the light of which it becomes clear 
that only vibrations of low frequencies may produce 
significant effects other than a simple renormaliza- 
tion of the pure electronic interactions. The model 
pertinent parameters are therefore U/t, V/t, lo jt 
and gjijj 



B. Computational details 

The calculations have been carried out using the 
infinite system density-matrix renormalization group 
(DMRG) methoclE3 with open boundary conditions. 

Since an infinite number of phononic quantum states 
lives on each site we have truncated the phonons basis set 
according to the previous section analysis, that is we kept 
only the two lowest vibronic states on each site (i.e. the 
two lowest eigenstates of the on-site Hamiltonian). This 
choice is physically reasonable since (i) we work at T = 



and therefore only the lowest vibronic states are expected 
to be involved, (ii) the molecules form well defined enti- 
ties that are only perturbatively modified by the presence 
of their neighbors. As already mentioned, when an elec- 
tron hopes between two nearest neighbors sites the hop- 
ping integral is renormalized by the Franck-Condon fac- 
tors, i.e. the overlap between the initial and final vibronic 
states of the sites, {{n,v,Sz)i ; (n', i/', S'z')i+i|_ff |(n — 
l,^M,Sz±l/2),■,{n' + l,^l',SzTlm^+l) =t{u\ti){v'W) 
(n and n' being the number of electrons on sites i and i+\ 
respectively, v and /i, v' and /i' the phonons states and 
Sz, Sz' the spin projection quantum numbers). Figure H^ 
shows the overlap between the vibronic ground state of 
a site supporting n electrons and the vibronic states of 
the same site supporting n ± 1 electrons. As can be seen, 
when g/uj is small the overlap, and therefore the Franck- 
Condon factor, decreases very quickly with the number 
of bosons, thus only the first few vibronic states arc in- 
volved and the truncation is totally pertinent. This fact 
is confirmed by exact calculations on small systems (4 
sites) where four phononic quantum numbers have been 
considered. For instance, the weight of the 3 and 4 bosons 
contributions in the wave functions is only 0.012 for 
g/to = 0.5, U/t = 4 and U /V = 4. When g/uj increases, 
the maximum of the Franck-Condon factors is rapidly 
displaced toward very large number of bosons. These 
vibronic states being strongly hindered by their large vi- 
brational energy, they have a small weight in the wave 
function and the system tends to localize. For instance 
for g/tu = 3, U/t = 4 and U/V = 4, the weight of the 3 
and 4 bosons contributions in the 4 sites system is only 
1.2 X 10~^. In the intermediate region the contribution 
of intermediately large phonons states (with occupations 
3,4,5) is not as negligible (0.3 for g/uj = 1.5, U/t = 4 
and U/V = 4), and the basis set truncation lowers the 
total hoping between nearest neighbors sites and thus in- 
creases the system localization. One can therefore expect 
that in the intermediate regime the true phases transi- 
tions will be displaced — compared to our results — 
toward larger values of the electron-phonon interaction. 
It is however clear that any truncated basis set will have 
a great deal of problem to accurately treat systems too 
close to a phases transition since the effective softening 
of the frequency near the transition translates into the 
implication of a quasi-infinite number of phonons states. 
From the above analysis, one can be quite confident in 
the quality of the numerical results and in particular in 
the different phases found in this work, provided that the 
exact position of the transitions is not seek at, in the in- 
termediate coupling region. In order to estimate more 
precisely the transition displacement due to the basis set 
truncation, we have run additional calculations using up 
to three bosons states per site occupation number (that 
is 12 on-site states) for a set of chosen electronic pa- 
rameters {U/t = 4, U/V = 4) and all values of the 
electron-phonon coupling constant. 

In order to characterize the phases diagram, we have 
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FIG. 1: Overlap between the vibronic ground state Vq of a 
site supporting n electrons and the a**^ vibronic state z/^ of 
a site supporting n ± 1 electrons . 



computed the charge and spin gaps, defined as usual : 

Ap = £;o(2iV, A^+1, 0) + £;o(2iV, A^-1, Q)-2Eq{2N, N, 0) 

and 

A,, = Eq{2N, N, 1) - Eq{2N, N, 0) 

where Eo{Ns,Ne,Sz) is the ground state (GS) en- 
ergy of a system of N^ electrons, Ng sites and spin 
projection Sz. In addition, we have computed the 
charge-charge, spin-spin correlation functions : ca{J) = 
{{Aq - (Aq)) {Aj - {Aj))) where A stands either for the 
number operator, n, or for the spin projection opera- 
tor, Sz, and the on-site singlet correlation function : 
csgij) = {Sgo'Sgj) where Sgi'^ = 4^4^ ^^d Sgi are 
the singlet creation and annihilation operators on a site. 

The properties calculation have been computed using 
255 states per renormalized block, whereas for the gaps 
calculations we have used a double extrapolation (i) on 
the systems size, and (ii) on the number, to, of states 
kept, using m = 100, 150 and 255. In all calculations 
we have computed systems of size up to 80 sites and 
extrapolated to the infinite chain. 

To have more information on the localized phases wave 
functions we have also computed the density matrices at 
the central sites and performed exact diagonalization of 
small systems. 

III. RESULT 

The present work explores the whole range of the on- 
site repulsion strength and of the e-mv coupling. The 
vibration frequency have been chosen to be uj/t = 0.2 
for the reasons already exposed in the preceding para- 
graphs. The phases diagrams have been computed for 
two values of the nearest-neighbor versus on-site repul- 
sions ratio V/U which are recognized to be generic for the 
ID organic conductorsEil, namely V = C//4 and V = U /2. 



Phase diagram 

V=U/4 £0=0.2 




OLL 


OLE 


A 4kp CDW 


A LRO 4kp CDVr 


T LRO 2kp CDVr 


+ BP 


*P 


— U,„/U=8.5 

- V,(N.=4) 



regime one has the polaronic (t4// > 0) and bi-polaronic 
{Ueff < 0) phases where the electrons are self-trapped 
(alone or by pairs) by the molecular geometry deforma- 
tions. In between these two regimes, that is for interme- 
diate e-mv coupling, one finds the two new phases. The 
first one is an insulating long-range ordered Akp CDW 
phase which develops at the extend of the metallic ikp 
CDW phase for Ueff > 0. The second one is an insulat- 
ing long-range ordered 2A;f CDW phase which develops 
for Ueff < at the expends of the localized bi-polaronic 
phase. It is noticeable the Ueff — line seems to remain 
a strict phases boundary. On the contrary the other fron- 
tiers have been shifted. For Ueff > the localized phases 
are enhanced and the delocalized ones reduced. On the 
contrary, for Ueff < the bi-polaronic phase is reduced. 
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FIG. 2: Phases diagram of the extended Hubbard-Holstein 
model for V = [7/4 and co/t = 0.2. 
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FIG. 3: Phases diagram of the extended Hubbard-Holstein 
model for V = U/2 and uj/t = 0.2. 



Figure H and |3| report the phases diagrams for V = U/A 
and V = U/2 as a function of g/oj and U /t. The two 
diagrams present the same general features with seven 
different phases. The major effect of the introduction of 
NN repulsion in the Hubbard-Holstein model is to sta- 
bilize two new phases in the e-mv intermediate coupling 
regime. From another point of view, the inclusion of 
the e-mv coupling in the extended Hubbard model has 
similar consequences as its inclusion in the pure Hub- 
bard model; that is : the apparition of polaronic and 
bi-polaronic phases in the strong coupling regime, the ap- 
parition of a Akp CDW phase in the intermediate regime 
for extremely low values of the on-site repulsion. 

To summarize, in the weak coupling regime one find 
both the Luttinger Liquid phase for Ueff > and the 
Luther Emery phase for Ueff < 0. In the strong coupling 



A. The Ueff > phases 

1. Luttinger liquid phase 

For small values of g/io, up to intermediate ones if the 
on-site repulsion U is not too large, one finds, both in 
the U/V = 4 and U/V — 2 cases, the expected LL phase. 
The computed charge and spin correlation functions ex- 
hibit power law behavior with dominant 2kp SDW fluc- 
tuations and sub-dominant CDW fluctuations. The spin 
and charge gaps extrapolate nicely to zero to numerical 
accuracy. 

Similarly _to what happens in the pure extended Hub- 
bard modeltj, the 2kF SDW fluctuations and the Akp 
CDW fluctuations are enhanced by increasing values of 
the NN repulsions. From the charge structure factor 
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FIG. 4: Kp as a function of g/uj for U = 0.2, 1, 2.5, 4, V = 
[7/4 (solid lines) and V — U/2 (dashed hnes). The symbols 
follow the same pattern as in the phases diagrams. 



Sp{q) we have computed the LL Kp parameter as 
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Figure 11 reports the Kp parameter as a function of both 
U/V and g/w. Once again the resuhs are a simple su- 
perposition of the Kp reduction effect due to the NN 
repulsions and the reduction effect due to the e-mv cou- 
pling. 



2. TheAkp CD W phase 

The LL phase is bordered, both for V — U /A and for 
V — f//2, by a metallic phase presenting dominating Akp 
charge fluctuations. This phase have very similar chat- 
acteristics as the Akp CDW phase found for V = 0li3, 
that is no charge neither spin gap (to numerical accu- 
racy), power law decreasing of the charge and spin cor- 
relation functions, dominant Akp CDW fluctuations and 
very small values of Kp compared to the purely electronic 
model. One has for instance, for U /t = 4 and V/t — 1, 
Kp = 0.28 when g/ui = 1.5 instead of Kp = 0.55 in 
the pure electronic cH model (see fig. H). One should 
however notice that while the Kp values remain always 
larger ptian the 1/4 minimal value predicted by the LL 
theoryll3 for metallic behavior, it can be as large as 0.48 
for U = 2.5, g/u = 1.25 and V = C//4, that is much 
above the 1/3 limiting value predicted, by the LL theory 
for dominant Akp CDW fluctuationstj. Despite the ab- 
sence of long range coulombic repulsions, this phase is in 
fact, in many ways very alike a Wigner crystal. 

Figures ^ and ^ show that the NN repulsions have a 
strongly destructive effect on this phase. Indeed, it re- 
duces strongly its domain of existence for increasing g/uj. 
Compared to the V = Q case a new insulating, long- 
range order (LRO) Akp CDW phase has taken a large 
part of the g/to parameters range of the metallic Akp 
CDW phase. For increasing V/U the metal to insulator 
phases transition (MIT) is repelled to smaller values of 
g/u}, squeezing the metallic Akp CDW phase toward the 
LL one. 



3. The LRO ^kp CDW phase 

As g/io increases, the system undergoes a metal to in- 
sulator phases transition and the Akp CDW fluctuation 
phase condensates into a long range ordered Akp CDW 
phase. 

In order to characterize this new phase, we have 
computed the staggered charge correlation functions 
{-lyC^iJ) where 

Cn{j) = (("i - n){n^+j - n)) 

and n = N^/Ns = 1/2 is the average charge per site. The 



associated order parameter is therefore 

X4,, = hm ^(-l)JC„(j) 

J 

In this gaped regime, one should be careful and clearly 
distinguish between the correlation functions C„(j) and 
the correlation functions of the observable fluctuations 
Cn{j). Indeed, while in delocalized phases the two do not 
differ, in gaped phases the correlation functions tends 
toward a non zero constant as the inter-site distance in- 
creases, while the fluctuatipas decrease quickly to zero at 
infinite inter-site distanceal3. 
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FIG. 5; (a) Order parameter X^tp and (b) charge gap as 
a function of g/uj for [7 = 4. The solid lines correspond to 
U/V — 2 and the long dashed lines to U/V = 4. The symbols 
follow the same pattern as in the phases diagrams. The short 
light dashed lines correspond respectively to Ap = V (that is 
V = 2 and l^ = 1) and Ap = U^ff. 

Figure o reports both the charge gap Ap and the order 
parameter X^kp as a function oi g/cu, for the two values of 
the NN repulsion. One sees immediately that the opening 
of the charge gap from the metallic Akp CDW phase is 
simultaneous with the formation of the long range order. 
Similarly the computed charge fluctuations correlation 
functions c„(j) go from a power law behavior as a func- 
tion of increasing inter-site distances to an exponential 
behavior. One should note that the spin channel remains 
ungaped and the corresponding fluctuations correlation 
functions cszij) decrease as a power law with increasing 
inter-site distances. 

At the metal to insulator phases transition, that is for 
intermediate values of g/uj (—1), one observes a smooth 
opening (exponential like) of both the charge gap and the 
order parameter ; opening that seems consistent with a 
Kosterlitz-Thouless transition. For large value of g/oj, 
the order parameter saturates and the system undergoes 
a self-trapping transition of the electrons toward a pola- 
ronic phase. It is noticeable that, while the MIT is very 



soft, the self-trapping transition is on the contrary rather 
sharp. 



4-. The small polaronic phase 



As expected from previous worksElllallZIO, in the strong 
couphng regime but stiU positive effective on-site repul- 
sion, Ueff — U ~ 2g'^ /oj, the system undergoes a tran- 
sition toward a polaronic phase where the electrons are 
self-trapped by the molecular distortions. This trapping 
is mediated by the Franck Condon factors, that strongly 
renormalize the hopping integrals between low energy 
vibronic states. One can remember that the hopping 
between the ground vibronic states is renormalizcd as 
t — > teff ~ t exp{—{g/uj)'^). The ground state of the sys- 
tem is therefore dominated by configurations such as : 

..10101010.. 

where the 1 stand for sites supporting one electron and 
the stand for empty sites. The validity of this pic- 
ture has been checked both on the GS wave function 
of small systems (4 and 8 sites, PBC, exact diagonal- 
ization) and the on-site density matrix in the DMRG 
calculations. On small systems the computed weight of 
those configurations is always larger 0.85 with, for in- 
stance, 0.900 ior U = A, V = 2, g/cu = 2 and 4 sites. On 
large systems, we have computed the central sites den- 
sity matrices and found that the probability of having 
double occupations is extremely small, with for instance, 
p(Ti) < 10-9 for C/ = 4, g/uj = 2.5 and aU values of 
U/V. Coherently, the GS energy per site is nearly inde- 
pendent of the system size and verify — up to at least 
4 significant numbers — on all the computed points, the 
formula —N^/Ns g'^/uj+uj/2. Such a GS is strongly quasi- 
degenerated due to the equivalence between the odd and 
even sites and the different spin configurations. The 
small splitting is due to the residual delocalization and 
therefore scales as te// = iexp{— (gi/w)^) (see figure ||). 
The spin channel remains ungaped. The main difference 
between the present phase and the phase found in the HH 
model stays in the charge channel. Indeed, the NN re- 
pulsion is responsible for the opening of a strong gap (see 
figure H) that did not exist in the V = case. In fact, 
the charge gap scales as the cost to add an extra elec- 
tron to the system. In the case of open systems the end 
sites being always occupied, Ap scales as min {Ueff,V) 
(according to whether the extra electron is located on an 
"empty" site or on an already "occupied" one) while in 
periodic systems is scales as imn {Ueff ,2V). The change 
of behavior in Ap can be clearly seen on figure ||, where, 
for instance, the gap for [/ = 4, U/V = 4 undergoes a 
saturation to Ap = V^ = 1 at the self-trapping transition 
(that occurs between g/uj = 2 and g/uj = 2.2) and then 
a strong decrease for g/uj > ^15/2 ~ 2.74 {Ueff = V), 
where is behaves as C/ — 2g'^/uj. One should notice that 
the full saturation of the order parameter occurs only 
after the second transition. 
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FIG. 6: (a) First excitation energy, Aio, as a function of 
g^ /lo^ for 4 sites systems, several value of V and for U = 41/ 
(circles), U = ZV (solid lines) and U = 2V (crosses). For 
the V — case we have chosen U — S, U = 4: and U — 2. 
(b) Phase boundary between the Akp LRO CDW phase and 
the polaronic phase as a function of V and teff. (c) Aio as a 
function of t^ff in a log-log curve for small values of teff (for 
clarity, we only show the curves for V = 2). 



In order to better study the position of the phases tran- 
sition between the Akp LRO CDW and the polaronic 
phase, we performed exact diagonalizations on periodic 
small systems (4 sites, 2 electrons). One should remem- 
ber that while the Akp LRO CDW phase has a non- 
degenerated GS, the polaronic phase presents a quasi- 
degenerated GS, the degeneracy hfting being of the or- 
der of magnitude of teff . We have therefore used the first 
excitation energy, Aio, as a criteria for the phases tran- 
sition. Figure ||(a) reports Aio as a function of g'^/u;'^ 
for different values of U and V. One first notice that the 
excitation energies depend in a negligible way on on-site 
repulsions. On the contrary it does strongly depend on 
the NN repulsion. Going from very large e-mv coupling 
to smaller values, the excitation energy first increases as a 
power law of tg// = texp(— g'^/w^) (see figure g(c)) in the 
polaronic phase then linearly as a function oi g^ /uj^ in the 
Akp CDW LRO phase. Decreasing g^ joj^ to even lower 
values this excitation energy should go through a max- 
imum and then decrease back to zero at the MIT. The 
location of the phases transition between the 4A:p LRO 
CDW and the polaronic phase have been evaluated as the 
point where the linearly extrapolated excitation energy 



of the 4kp LRO CDW phase crosses the zero axis. Fig- 
ure ||(b) reports the phase boundary as a function of t^ff 
and V. One sees immediately that it foUows a perfectly 
linear curve that can be fitted as T^ = 129.36ie// — 0.93i, 
in these coordinates. This curve have been reported on 
the phases diagrams (figures y and y) as the Vc curves. 
One sees immediately that the phases transition position 
is very weakly dependent of the system size (as expected 
from such localized systems) and that the small systems 
estimations work pretty well for the infinite systems. 



B. The Ueff < phases 

1. The Luther-Emery phase 

For negative values of Ueff and small values of g/to, 
we found a metallic phase for which all spin, charge and 
on-site singlet fluctuations correlation functions decrease 
with the inter-site distances, as power law. All three cor- 
relation functions exhibit dominating 2kp components. 
In all computed cases, the 2fci? CDW fluctuations have 
the largest amplitudes. The main effect of the NN repul- 
sion is to increase the amplitude of the CDW fluctuations 
and to strongly decrease the amplitude of the on-site sin- 
glet fluctuations. The charge gaps clearly extrapolate to 
zero, whereas we found a very small gap in the spin chan- 
nel (A^ ~ 0.002 - 0.003). This fact is not incompatible 
with the behavior of the spin-spin correlation function, 
since very small gaps means very large correlation lengths 
of the order of magnitude of A~^. The expected expo- 
nential behavior of the spin correlation functions should 
therefore take place at inter-site distances larger than 
the computed chain lengths. One can easily recognize 
in this phase a weakly gaped Luther-Emery phase. The 
values of Kp extracted from the charge structure factors 
are strongly reduced compare to the values of the purely 
electronic model, in a similar way as what has been found 
in the V = case. It should be noted that Kp always 
remains lower than 1 in agreement with dominant CDW 
fluctuations. 



2. The LRO 2kF phase 

When g/iu increases toward the intermediate regime 
{g/uj > 1.5), the system undergoes a MIT toward an in- 
sulating phase presenting a 2kp long range order. One 
should notice that this phase is induced by the NN repul- 
sions and does not exist in the Hubbard Holstein model. 
In comparison to the V ~ case, the 2kp LRO phase 
develops at the expense of the bi-polaronic phase . It is 
interesting to point out that for positive Ueff the devel- 
opment of the 4:kp LRO phase, induced by the NN repul- 
sions, have a tendency to localize the electronic structure, 
while for negative Ueff , the development of the Akp LRO 
phase corresponds to a tendency toward a less localiza- 
tion. The amplitude of the charge correlation functions 



Cn{j) extrapolate at infinite inter-site distances toward 
finite values (for instance 0.06 for C7 = 1, y = 0.5 and 
g/u! = 2). The order parameter has been defined in the 
usual way 
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and is reported on figure m. One can see that, as in 
the Akp LRO phase, the order parameter increases very 
slowly at the MIT as in a Kosterlitz-Thouless transition. 
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FIG. 7: Order parameter X2fcp as a function of g/uj for U = 1. 
The solid lines correspond to U/V = 2 and the long dashed 
lines to U/V = 4. The symbols follow the same pattern as in 
the phases diagrams. 

Both spin and charge channels are gaped. It is notice- 
able that the gaps values are always of the same order of 
magnitude : Ap = A^ = 1.04 for U = 0.2, V == 0.1 and 
g/uj = 2, Ap = 0.35, A^ = 0.33 for C/ = 1, V" = 0.25 
and g/uj — 2. Coherently, the fluctuations correlation 
functions (spin, charge and on-site singlet) decrease ex- 
ponentially with the inter-site distances (see figure pi). 

One should note that, for large |f/e//|, the on-site sin- 
glet fluctuations correlations are dominant, whereas, for 
weak |C/e//|, the charge fluctuations correlations domi- 
nate. On the contrary the increase of V, increases the 
CDW at the expense of the singlet fluctuations. 



3. bi-polaronic phase 

For large values of g/uj the system goes in a bi- 
polaronic phase where the attractive nature of the ef- 
fective on-site interaction strongly couples the electrons 
in pairs. This pairing is associated with a strong local- 
ization and a self-trapping of the pairs, due to the rescal- 
ing of the hopping integrals between low energy vibronic 
states by Franck Condon factors. The remaining delo- 
calization processes are very small and the ground state 
wave functions have dominant configurations of the type 
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FIG. 8; Charge (circles) and singlet (crosses) fluctuations cor- 
relation functions for U/V = 2, g/ui = 2 and different values 
of U^ff. U^ff = 0.6 (a), U^ff = 1.1 (b) and U^ff = 1.4 (c). 



states instead of eight. Since these calculations are very 
expansive we have only performed them for a fixed set 
of electronic parameters U/t = 4 and V/U = 1/4. Vary- 
ing the electron-phonons coupling parameter therefore 
leads us to go from the Luttinger Liquid phase up to the 
polaronic phase. Figure |9| reports the order parameter 
Xikp for the two and three phonons states. As expected 
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where 2 stands for a double occupancy of a site and for 
an empty site. The probability of having a lonely elec- 
tron on a site is extremely small with for instance values 
smaller than lO""^ for U/t = 1, V/t = 0.25 and g/uj = 2.5. 
The energy per site is nearly constant and in excellent 
agreement with the formula l/AUeff — l/2g^/uj + 1/2 w 
(the difference between the computed values and the for- 
mula being of the oder of 10"^ for U/t = 1, V/t = 0.25 
and g/oj = 2.5). Let us notice that the GS is strongly 
quasi-degenerated due to the equivalence between the 
four different phases of the CDW and to the absence 
of second neighbors repulsion terms. As in the polaronic 
phase the degeneracy splitting scales as ie//- Finally the 
system is strongly gaped both in the charge and spin 
channels. Indeed, either to extract an electron from the 
system or to build a triplet state, one needs to break an 
electron pair. Such a mechanism cost the energy Ueff 
and therefore both gaps scale as it. 



C. The basis set truncation effect 

As we already mentioned the phononic basis set trun- 
cation to two vibronic states per occupation number 
should induce a bias toward excessive localization in the 
intermediate regime. This is precisely the range of cou- 
pling parameters where the new phases appears and, as 
we will see later, the coupling regime which is the most 
interesting for the physics of real systems. In order to 
quantify the effects of the basis set truncation we have 
performed, in additi on t o the small systems calculations 
presented in section II B , infinite system DMRG calcula- 



tions with three vibronic states kept for each site occu- 
pation and spin. Each site is then described by twelve 



FIG. 9: X4kp order parameter for two (upper curve) and 
three (lower curve) phonons states kept per electronic state 
on each site. The phases symbols follow the same pattern as 
in the phases diagrams. 



the increase of the basis set results in a displacement of 
the phases transitions toward larger value of the coupling 
constant, that is a larger delocalization of the system for 
a given value of g/tu. The Luttinger Liquid phase is only 
slightly enlarged, in agreement with the weak participa- 
tion of excited vibronic states in the small systems wave 
functions. The phase which is the most favored by the 
basis set increase is the metallic ikp CDW phase. Indeed 
this phase is substantially extended toward larger values 
of g/u! while the insulating LRO Akp phase is essentially 
translated by 0.75g/uj. 

In conclusion the basis set truncation does mot change 
the structure of the phases diagram, its main effect being 
of reducing the range of the metallic ikp CDW phase, 
essentially on the side of the larger values of the electron- 
phonons coupling parameter, and shifting the LRO 4/ci^ 
phase toward smaller values of g/uj. All these results 
are in complete agreement with the small system calcu- 
lations. 



IV. RELEVANCE TO THE ORGANIC 
CONDUCTORS PHYSICS 



In this section we are going to discuss the perti- 
nence of the above calculations for the low energy 
physics of organic conductors and more specifically for 



the {DX - DCNQT)^ M family and the Bechgaard salts 
family. Indeed, as already mentioned in the introduction, 
the building molecules of these systems present a certain 
number of similar characteristics that are crucial for the 
relevance of the molecular vibrations to the electronic 
structure. The DX - DCNQI as well as the TMTTF or 
TMTSF molecules are planar, strongly conjugated and 
their skeleton is based on organic cycles ; the quinone 
cycle for the DX — DCNQI and the two pentagonal 
cycles of the fulvalene for the TMTTF and TMTSF. 
These geometrical properties strongly favor the existence 
of low frequency vibrational modes, namely the angular 
distortion of the cycles, the bond length remaining un- 
touched. Among these, the Ag modes couple to the elec- 
tronic structure in a relatively strong manner, leading 
to e-mv coupling constants in the intermediate range (as 
can be seen in Raman experimentgj or in the gepmetry 
relaxation of the molecules with their ionizationlij). In 
addition, one should notice that due to their conjugated 
character, the tt system of the considered molecules are 
strongly delocalized on their skeleton and thus strongly 
polarizable. The consequence is that the second neigh- 
bor coulombic interactions should be strongly screened 
by the "metallic plate" of the in-between molecule. One 
can therefore reasonably assume that the on-site and first 
neighbor repulsion terms are the only pertinent ones for 
the physics of these organic ID systems. 

Recently charge ordered phases presenting similar 
characteristics with the Akp LRO CDW state have been 
observed both in the {DX - DCNQI)^ M family and in 
the Bechgaard salts family. 

The most characteristic compound is certainly the 
{DI — DCNQI)^ Ag compound. Indeed, this system 
which is strongly one dimensional undergoes a MiT 
phases transition at 22Cljr toward a 'ikp CDW stateE3. 
Both NMRB and X-rayB show that the insulating state 
exhibit a on-site charge disproportionation between two 
adjacent molecules. This charge order saturates at a 3 : 1 
ratio below lAQK and is associated with a molecular 
geometry deformation due to the ionicity modification. 
In addition, the spin channel remains ungapedEll. This 
compound, which is usually considered as strictly non- 
dimeuzed (despite a recent doubt raised by Meneghetti 
et ao) seems to be the perfect example of an e-mv driven 
Akp LRO CDW state as the one discussed in this paper. 
Most authors assume that this LRO state is driven by 
a Wigner crystallization due to strong long range elec- 
tronic repulsions. We have seen in the previous sections, 
that the e-mv driven Akp CDW phases present very sim- 
ilar characteristics as the Wigner crystals, without the 
need of very strong correlation effects and without long 
range repulsions. In the light of the previous consid- 
erations on the strongly polarizable electronic structure 
of the quinone cycles, on the usually assumed values 
of the electronic correlation strength (in the intermedi- 
ate rapge rather than the strong range for the DCNQI 
familyc^), and on the e-mv coupling characteristics, it 
seems to us more plausible that the considered Akp CDW 



state is driven by the e-mv than by the usual unscreened 
strong coulombic repulsions. At very low temperatures 
the [DI — DCNQI)^ Ag undergoes a second phases tran- 
sition toward a spin ordered state where the Akp LRO 
CDW is associated with a 2kp anti- ferromagnetic LRO. 
These results are in agreement with the 2kp SDW fluctu- 
ations exhibited in the Akp LRO CDW phase of the eHH 
model, fluctuations that could easily be pinned at low 
temperatures by impurities or inter-chains interactions. 

Another family for which the preceding analyses are 
relevant about the e-mv influence on the electronic struc- 
ture, is the Bechgaard salts family. This fact should 
be put in the light of the newly discovered Akp charge 
ordered phase. This phase have been observed on sev- 
eral systems such_as in the {TMTTF)^PFq or the 
{TMTTF)^ AsFq y. The charge ordered (CO) phase 
appears when the temperature is lowered from the metal- 
lic Luttinger Liquid phase through a soft cross over. It 
is noticeable that transport measurements observe Akp 
CO fluctuations in the normal LL phase, indicating that 
the interactions driving the transition are relevant far 
insid e th e normal phase and over a large range of pres- 
suresE§Ell. X-ray diffraction measurements do not see 
any n-merization associated with this phases transition, 
excluding a Peierls mechanism. NMR measurements on 
the carbon atoms of the central double bond of the ful- 
valene (on which the Highest Occupied Molecular Or- 
bital, responsible for the low energy physics of these 
systems, have large coefficients) exhibit a charge dispro- 
portionation on NN molecules. Magnetic susceptibility 
measurements are transparent to this phases transition 
and the spin channel remain ungaped. It is clear that 
all these experimental results are in full agreement with 
the characteristics found in the present work for the Akp 
LRO phases transition. Despite the fact that the sys- 
tem dimerization has not been taken into account in the 
present work, the e-mv coupling appears as a strong can- 
didate for the CO driving mechanism. It is however clear 
that the dimerization degree of freedom should be taken 
into account, in addition to the correlation degrees of 
freedom (both on-site and inter-sites) and the e-mv de- 
grees of freedom treated in the present work, for a com- 
plete description of the Bechgaard salts. 

Finally, one can point out that the mechanism under- 
lying a MIT transition toward an on-site charge modu- 
lation, driven by the e-mv, does not imply any lattice 
distortion apart from small modifications in the molec- 
ular geometries, in particular in the angles of the cycles 
(the fulvalene cycles in the Bechgaard salts). Indeed such 
geometrical relaxations can be expected to be a conse- 
quence of the charge disproportionation. These transi- 
tions should therefore be g = transitions. 



V. CONCLUSION 

The present paper studies the influence of the electron- 
molecular vibration coupling on the electronic structure 
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of correlated ID quarter- filled chains. The model chosen 
is the extended-Hubbard Holstein model, that is the sim- 
plest one containing both electron correlation effects rele- 
vant for the physics of molecular crystals such the organic 
conductors, and e-mv coupling whose effects can be ex- 
pected to play a crucial role. We have found that for low 
phonons frequencies the electronic structure is strongly 
affected by the presence of molecular vibrations. One 
of the most striking results being the existence of ikp 
CDW phases, one metallic and one insulating charge or- 
dered phase, for small values of the correlation strength 
(as small as U/t = 2), small values of the nearest neigh- 
bor repulsion and in the absence of long range coulombic 
repulsion. A study of these phases shows however that 
they have similar characteristics as a Wigner crystal. In 
this light a new interpretation — based on the e-mv cou- 
pling — of the origin of the Akp charge ordered phase 
in the {DI — DCNQI)2Ag and of the recently discovered 
CO phase in the {T MTT F)2X family has been proposed. 
Such an interpretation of the apparition of the Akp CO 
phase present the advantage over the usual purely elec- 



tronic interpretation to be in agreement with the usual 
values of the correlation strength for these systems and 
the strong screening of the long range bi-electronic repul- 
sions that can be expected from the tt conjugated char- 
acter of the building molecules. To conclude we would 
like to point out that for a complete description of the 
Bechgaard's salts physics, one should treat in addition to 
the intra-molecular phonons modes, the inter-molecular 
modes responsible for the known dimerization in these 
compounds. In|4dew of the recent work of Campbell, Clay 
and Mazumdaic3 on the adiabatic dimerized extended- 
Hubbard Holstein model that forecasts the existence of 
mixed phases such as the Bond Charge Density Wave 
and the Spin-Peierls ^kp CDW phase, it would be of in- 
terest to conduct the same type of study as the present 
one on the dimerized extended-Hubbard Holstein model. 
Indeed, such a model would include all degrees of free- 
dom important for ID organic conductors as well as the 
phonons quantum fluctuations that have been proved jta 
be crucial for a correct description of the e-mv couplingE^I 
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